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; 1 INTRODUCTION 

> ■ 

, The primordial cosmic gas is composed by Hydrogen (about 75 per 
» I ■ cent by mass), Helium (~ 24 per cent), and traces of other light el- 

' ements, like Deuterium, Helium-3, Lithium and Berillium. When 
this gas collapses into the dark matter halos typical of galaxy clus- 
ters (> 1O 14 M0), it undergoes shocks and adiabatic compression, 
reaching densities of about 10~ 3 particles cm -3 and temperatures 
of the order of 10 8 K. For these physical conditions, the plasma 
is optically thin and its X-ray emission is dominated by thermal 
bremsstrahlung processes. Its enrichment by metals (i.e. elements 
with atomic number larger than 5) proceeds by the releases to the 
medium of the products of the star formation activities that take 
place in the member galaxies. 

An effective way to estimate the metal abundance by num- 
ber relative to Hydrogen is to measure the equivalent width of the 
emission lines above the X-ray thermal continuum. Nowadays, the 
Iron abundance is routinely determined in nearby systems thanks, 
in particular, to the prominence of the K-shell (i.e. the lower level 
of the electron transition is at the principal quantum number n — 1) 
Iron line emission at rest-frame energies of 6.6-7.0 keV (Fe XXV 
and Fe XXVI). Furthermore, the few X-ray galaxy clusters known 
at z > 1 have shown well detected Fe line, given sufficiently long 
(> 200 ksec) Chandra and XMM-Newton exposures (Rosati et al. 
2004, Hashimoto et al. 2004, Tozzi et al. 2003). It is more difficult 
is to assess the abundance of the other prominent metals that should 
appear in an X-ray spectrum at energies (observer rest frame) be- 
tween ~ 0.5 and 10 keV, such as Oxygen (O VIII) at (cluster rest 



ABSTRACT 

We use models of the rates of Type la supernovae (SNe la) and core-collapse supernovae, built 
in such a way that both are consistent with recent observational constraints at z < 1.6 and 
can reproduce the measured cosmic star formation rate, to recover the history of the metals 
accumulation in the intra-cluster medium. We show that these SN rates, in unit of SN number 
per comoving volume and rest-frame year, provide on average a total amount of Iron that is 
marginally consistent with the value measured in galaxy clusters in the redshift range 0—1, 
and a relative evolution with redshift that is in agreement with the observational constraints 
up to z Ri 1.2. Moreover, we verify that the predicted metals to Iron ratios reproduce the 
measurements obtained in nearby clusters through X-ray analysis, implying that (1) about 
half of the Iron mass and > 75 per cent of the Nickel mass observed locally are produced by 
SN la ejecta, (2) the SN la contribution to the metal budget decreases steeply with redshift 
and by z « 1 is already less than half of the local amount and (3) a transition in the abundance 
ratios relative to the Iron is present between redshifts ~ 0.5 and 1.4, with core-collapse SN 
products becoming dominant at higher redshifts. 
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frame) 0.65 keV, Silicon (Si XIV) at 2.0 keV, Sulfur (S XVI) at 2.6 
keV, and Nickel (Ni) at 7.8 keV. 

A direct application of the detection of emission lines from 
these highly ionized elements is the possibility to describe the way 
the intracluster medium (ICM) is enriched of metals, whether by 
products from explosions of supernovae (SN) type la, mainly rich 
in Fe and Ni, or by outputs of collapsed massive stars (Core Col- 
lapsed SN), with relative higher contributions of q— elements like 
O, Si, S. However, the ability to resolve and measure elemental 
abundances other than Iron, which was observed in the early days 
of the X-ray analysis of galaxy clusters (e.g. Mitchell et al. 1976, 
Serlemitsos et al. 1977) has required a significant improvement in 
the spectral resolution and sufficiently large telescope effective ar- 
eas that have been reached initially with the X-ray satellite ASCA 
(e.g. Mushotzky et al. 1996, Fukazawa et al. 1998, Finoguenov, 
David & Ponman 2000) and then with the observatories of the new 
generation, such as Chandra (e.g. Ettori et al. 2002, Sanders et al. 
2004) and XMM-Newton (e.g. Bohringer et al. 2001, Gastaldello & 
Molendi 2002, Finoguenov et al. 2002). 

In the present work, we reverse the problem and infer from ob- 
served and modeled SN rates the total and relative amount of metals 
that should be present in the ICM both locally and at high redshifts. 
Through our phenomenological approach, we adopt the models of 
SN rates as a function of redshift that reproduce well both the very 
recent observational determinations of SN rates at z > 0.3 (Dahlen 
et al. 2004, Cappellaro et al. 2005) and the measurements of the 
star formation rate derived from UV-luminosity densities and IR 
data sets. We then compare the products of the enrichment process 
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Table 1. Adopted values for the average atomic weight (W), solar abun- 
dance by number with respect to Hydrogen (A, from Anders & Grevesse 
1989), total synthesized isotopic mass per SN event (rrij a and mcc from 
Table 1 in Nomoto et al. 1997; nucleosynthesis products of SN la are those 
from the deflagration model W7, whereas SN CC yields are integrated over 
the mass range 10 — 50AfQ with a Salpeter IMF) and corresponding abun- 
dance ratios by numbers with respect to the Iron, with each abundance nor- 
malized to the solar photospheric value A (Yi = m,i/m Fe X W Fe /Wi X 
A Pe /Ai). 
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Fe 


55.845 


4.68e-5 


0.743 




0.091 




O 


15.999 


8.51e-4 


0.143 


0.037 


1.805 


3.818 


Si 


28.086 


3.55e-5 


0.153 


0.538 


0.122 


3.526 


S 


32.065 


1.62e-5 


0.086 


0.585 


0.041 


2.284 


Ni 


58.693 


1.78e-6 


0.141 


4.758 


0.006 


1.647 



to the constraints obtained in the X-ray band for galaxy clusters 
observed up to z < 1.2. We note that previous work on the pro- 
duction of Iron in galaxy clusters by, e.g., Matteucci & Vettolani 
(1988), Arnaud et al. (1992), and Renzini et al. (1993) made use of 
SN rates that were known, with large uncertainties, only at z ~ 0. 
However, they concluded that either the SN la rates were higher 
by a factor of 10 in the past or the production of Iron in clusters 
by core-collapsed supernovae (SNe CC) would require a very flat 
Initial Mass Function (IMF). An alternative scenario is that SN CC 
produce 1/4 of the Iron and the remnants 3/4 is released from SN 
la, which would also be in accordance with the standard chemical 
model for the galactic halo and disk, allowing an IMF in elliptical 
galaxies similar to the one in the solar neighborhood. 

The main purpose of this work is to infer the mass of the most 
relevant and X-ray detectable elements present in the ICM and their 
relative abundance as function of redshift by using reliable (i.e. in 
agreement with several different observational constraints) models 
of SN explosion rates that include a dependence on the cosmologi- 
cal time. The fact that these rates are in unit of number per comov- 
ing volume and rest-frame time allows us to use them also over cos- 
mological distances typical of high-redshift systems. The method 
presented here is a first attempt to provide well-motivated and ana- 
lytic predictions of the expected metal abundances at z > 0.5 and 
of the relative role played by SNe la and SNe CC in polluting the 
ICM over time. 



2 OBSERVATIONAL CONSTRAINTS ON THE METAL 
CONTENT OF GALAXY CLUSTERS 

The observed Iron mass is obtained as 

M Fo , obs = 4nA Fe W Fe / Z Fe (r)p H (r)r 2 dr, (1) 
Jo 

where Z Fa (r) is the Iron abundance relative to the solar value 
A Fe measured in the ICM as function of radius, pu(r) = 
Pgas(r)/(2-2/i) is the Hydrogen density, p is the mean molecu- 
lar weight and has a value of 0.60 for a typical fully ionized cluster 
plasma with a number density n gas = nu + n c = 2.2nu. Our es- 
timates refer to a cosmology with parameters (Ho, f2 m , ^a) equal 
to (70,0.3,0.7). 

De Grandi et al. (2004) have correlated the total amount of 
M Fe observed with the X-ray satellite BeppoSAX in 22 nearby 



hot (fcT > 2 keV) galaxy clusters with the gas tempera- 
ture. For the assumed cosmology, we obtain a best-fitting re- 
sult of log(A/ Pc /lO lo Af ) = 0.33(±0.07) + 1.62(±0.53) x 
log(fcT/5keV) (scatter of 0.28). This fit considers quantities es- 
timated at the radius R500, where the cluster overdensity is 500 
times the critical density for an Einstein-de Sitter universe (i.e. an 
overdensity of 285 at 2 = for the cosmology assumed here) 
and corresponds to about 0.6 x R v i T . The overall results pre- 
sented in this paper do not depend upon the cluster gas temper- 
ature, although a value of 5 keV is hereafter adopted as refer- 
ence. Therefore, an Iron mass of (2.1 ± 0.3) x 1O 1O M is as- 
sociated with a 5-keV cluster within R500 at a median redshift 
of 0.05, with a scatter in the range of (1.1,3.9) x 10 10 M Q . At 
redshift between 0.3 and 1.3, for a sample of 19 objects observed 
with Chandra and presented in Tozzi et al. (2003) and Ettori et al. 
(2004a), that have a detection significant at the level of 2a of the 
Iron line emission and for which a single emission-weighted tem- 
perature and the radial profile of the gas density were measured, 
we obtain a best-fitting log(M Fe /lO lo M ) = 0.66(±0.06) + 
1.47(±0.48) x log(fcT/5keV) (scatter of 0.23) and infer a to- 
tal Iron mass for a 5 keV cluster at the median z — 0.63 of 
M Fo (< R 50 o) = (2.7 ± 0.4) x 10 10 M Q , with a scatter limited 
within the values of (1.6,4.6) x 1O 1O M0. These values are plot- 
ted in Fig. Q We assume that most (if not all) of the mass of the 
metals is located within -R500, even though this radius encloses 
just about 13 per cent of the virial volume (R500 ~ 0.5R v i T ). 
A further contribution from the outer cluster regions is however 
marginal: from numerical simulation (e.g. Ettori et al. 2004b) and 
extrapolated data (e.g. Arnaud, Pointecouteau & Pratt 2005), one 
can estimate that M5oo/Af v i r « 0.6 and, by assuming that in the 
outskirts (1) the metallicity is < 0.1 times the solar value, Aq, 
whereas its mean value in the central regions is about 0.3Aq, and 
(2) the gas follows the dark matter distribution, one can evaluate 
thatM Fo (i?500 < r < R vil )/M Fc (r < R 50 o) < 0.2. 

Apart from Fe, the most prominent lines detectable in X-ray 
spectra, like O, Si, S and Ni, were investigated initially with data 
from ASCA (an analysis of the average abundance of all the galaxy 
clusters in the ASCA archive is presented in Baumgartner et al. 
2005) and more recently with the larger effective area and improved 
spectral and spatial capabilities of XMM-Newton (e.g. Tamura et al. 
2004 present statistics of abundance ratios, also resolved spatially, 
for a sample of 19 nearby galaxy clusters). In Fig.|2| we summarize 
these results by plotting the relative number abundance ratios. The 
Iron abundance is about 0.3 A©, although the ASCA measurements 
have a very narrow distribution around the mean of 0.27 (±0.01 
at 90 per cent level of confidence) that departs significantly from 
the spatially resolved XMM-Newton estimates (only in the outer ra- 
dius, between 200 and 500 h^ 1 kpc, where the Fe/H decreases by 
about a factor of 2 with respect to the central estimate, it becomes 
consistent with the ASCA mean with a value of 0.32 ± 0.08). More- 
over, the O/Fe ratio is about solar, whereas the Si/Fe and S/Fe ra- 
tios are roughly super-solar and sub-solar, respectively. The average 
Ni/Fe ratio in the ASCA sample is about 3.4Aq, well in agreement 
with the original ASCA determination in the data of the Perseus 
cluster by Dupke & Arnaud (2001) and consistent with the XMM- 
Newton constraints in Gastaldello & Molendi (2004) and with the 
BeppoSAX results for a sample of 22 objects (Fig. 6 in De Grandi 
& Molendi 2002). On the other hand, it is worth mentioning that 
the measurements of the Nickel abundance becomes reliable only 
in hot (kT > 4 keV) systems where the excitation of (mainly) K- 
shell lines makes them detectable in the spectrum around 8 keV, 
in a region less contaminated from blends with Iron but highly af- 
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Figure 1. Accumulation history of the Iron mass from zp = 10 to z = 
0, as obtained from equations (2) and (3). The solid line shows the sum 
of SN la (dashed line) and SN CC (dotted line) contribution. Horizontal 
lines indicate the redshift region where the SN rates are actually measured 
(Dahlen et al. 2004). The accumulation history is here compared to the total 
Iron mass measured in samples of local (De Grandi et al. 2004) and high— z 
(Tozzi et al. 2003) galaxy clusters for a typical object at 5 keV. Error bars 
on the mean observed value and scatter (shaded region) around it are shown 
(see text in Section 2). The dot-dashed line indicates the expected Iron mass 
for a typical 5 keV cluster with M v j r = 10 15 Mq, a gas mass fraction of 
10 per cent and Iron abundance of 0.3^4.0. 



fected from a proper background subtraction. These difficulties put 
the actual limitations on solid Ni measurements for a large data set. 

Given the differences between the two data sets, and the fact 
that the XMM-Newton estimates are spatially resolved and thus sen- 
sitive to the presence of any gradient that instead is washed out in 
the global measurement shown here from ASCA, they are consid- 
ered independently in the following analysis. 



3 ACCUMULATING THE METALS 

The total predicted Iron mass is obtained as 



M Fe ,Ia + M Fe ,cc 

mFe.Ia x ria(di) x dt x dV 



dt,dV 

E 

dt,dV 



+ 



niFe,cc x r cc(dt) x dt x dV 



(2) 



where we assume mp e ^. 
ble 1 (from Nomoto et 

at given redshift equals to dt(z0, zl) = Hq j" (1 + 



and rriFe,cc as quoted in Ta- 
al. 1997) and an interval time 

zl , 



1 [n m (i + z ) 3 + (i - 



fin 



q a ){i + z) 2 + n A 



-0.5 



dz. 



As cluster volume, we define the volume corresponding to 
the spherical region that encompasses the cosmic back- 
ground density, p h = 3H$/(8nG) x Cl m (k 4x 1O 1O M 



Figure 2. Predictions (big dots) of the local metal yields compared with the 
X-ray measurements in Tamura et al. (2004; case for a medium-T cluster 
in Table 6 with estimates in the three radial bins — 50, 50 — 200 and 
200 - 500ft" 1 kpc; diamonds) and Baumgartner et al. (2005; 4.5 keV bin 
in Table 4; triangles). Error bars at lcr level are indicated. The abundance 
measurements are relative to the Anders & Grevesse (1989) solar photo- 
spheric values quoted in Table 1 . 



Mpc 3 for the assumed cosmology), with a cumulative mass of 
M vir » 6.8(fcT/5keV) 3/2 x lO 14 Af , as inferred from the best-fit 
results on the observed properties of X-ray clusters with kT > 3.5 
keV in Arnaud et al. (2005) and adopting M 500 /M vil « 0.6: 



dV = Afvir/Pb 



24, 514 



,3/2 



Mpc 3 . It is 



10 1B M o V5keV7 

worth noticing that the A/ V i r — T relation for simulated galaxy 
clusters (e.g. Ettori et al. 2004b), with an emission-weighted 
temperature greater than 3.5 keV, provides a typical virial mass 
(and corresponding volume) that is larger than the above value 
by about 20 per cent. We discuss in Section 4 the effect of this 
discrepancy on the overall results. The volume measured at z = 
is considered as the dimension of the region involved in the 
formation of a typical galaxy cluster and is maintained constant in 
redshift. 

The supernovae rates, i.e. the number of SN per rest-frame 
time and comoving volume, come from models that fit well the 
observed values and reproduce properly the Star Formation Rate 
(SFR) as widely discussed in Strolger et al. (2004) and Dahlen et 
al. (2004). In summary, once a functional form is defined to de- 
scribe the evolution with time of the SFR (see equation 5 in Strol- 
ger et al. 2004 with parameters for the extinction-corrected model, 
called M 1 ; the uncorrected M2 model just reduces the total amount 
of metals released by a factor of ~ 4), the corresponding SN rates 
can be defined as 



ri»(t) = v I SFR(t') x 4>(t - t') dt' 



rcc(i) 



kxh 2 x SFR(t) 



(3) 



where the normalization v indicates how many SNe la explode per 
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Figure 3. SN la rates and number ratio of SNe CC to la as a func- 
tion of redshift for the delay time distribution functions considered in the 
present work. (Solid line) A "narrow" Gaussian with r = 4 Gyr and 
ut d = 0.2t, which is the one adopted here; (dashed line) a "wide" Gaus- 
sian with r = 4 Gyr and at . = 0.5r; (dotted line) a exponential form, 
<fi(td) oc exp(— t^/r), with r = 5 Gyr. A similar figure is presented and 
discussed in Dahlen et al. (2004). 
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Figure 4. Metal abundance ratios as a function of redshift. Horizontal 
lines indicate the redshift region where the SN rates are actually measured 
(Dahlen et al. 2004). 



unit formed stellar mass and is function of the delay time distribu- 
tion function <j>(td) that represents the relative number exploded at 
a time td since a single burst of star formation; tF is the time of 
formation of the first stars and corresponds to zf = 10; the nor- 
malization k is the number of CC progenitors per unit of formed 
stellar mass; h is the Hubble constant in unit of 100 km s -1 Mpc -1 . 
In the following analysis, we adopt a "narrow" Gaussian form for 
with a delay time r = 4 Gyr and a td = 0.2r that better 
fit both the local (z < 0.1) rates and the GOODS data (Strongler 
et al. 2004, Dahlen et al. 2004) and fix v = 0.0010. We assume 
k — 0.0069MQ 1 as appropriated for CC progenitors masses in 
the range 8M < M < 50M© with a Salpeter (1955) Initial Mass 
Function (IMF). We show in Fig.[3]the relative rates of SNe CC and 
la as function of redshift: this ratio ranges between 2 and 3 up to 
z ~ 1 and increases steeply at higher redshifts, with the "narrow" 
Gaussian function that provides a SN la rate at z ~ 2 lower by an 
order of magnitude than the other delay time distribution functions. 

Regarding the adopted delay time, we mention here that Gal- 
Yam & Maoz (2004) and Maoz & Gal- Yam (2004) support the ar- 
gument that some current data on the cosmic (field) star formation 
history require for the observed rate of SNe la a delay time larger 
than 3 Gyr. Combining this with the low observed rate of cluster 
SN la at z < 1, they conclude that the measured amount of Iron in 
clusters cannot be produced mainly by SNe la explosions but, more 
probably, must be produced via outputs of SNe CC originating from 
a top-heavy IMF. The reader is refereed to the above-mentioned 
work for a detailed discussion on the effects of the assumed IMF 
on the SN rates. Our approach here is to use the best phenomeno- 
logical models of cosmic star formation history to infer the corre- 
sponding cluster metal accumulation history and to compare some 
predictions with the observational constraints available. 

To produce the cluster metal accumulation history, we have 
to take into account the fraction of the total produced Iron mass 
that remains locked up in the stars. The fraction of the total Iron 
that is then released to the ICM is estimated by considering that 

M s tar. We adopt Zicm ~ O.SAq 
(see Fig.|2j, M s tar ~ 0.13Micm (from equations 1 and 10 in Lin, 
Mohr & Stanford 2003 for a 5 keV cluster) and Z Bta r = 1.3&4 
(from interpolation of the results in Lin, Mohr & Stanford 2003; see 
also Renzini 2003, Portinari et al. 2004). Therefore, Mf c ,icm = 
Mp e ,tot/(l + fz), where f z = M st8r /MicM x Z stSjI /Zi C M » 
0.59. Considering its weak dependence upon time (see, e.g., Ta- 
ble 2 and Fig. 3 in Portinari et al. 2004), we neglect the evolution 
of the locked-up fraction. 

In Fig. □ we plot the results of the expected Iron mass com- 
pared to the observational results discussed in Section 2. The pre- 
dicted amount is well in agreement with the observed Iron mass 
in local clusters. At a median redshift of 0.05, the De Grandi et 
al. sample has a mean Mp e associated to a 5 keV object that is be- 
tween 30 (when the M v i r — T relation from simulations is adopted) 
and 60 (with the observed A/ V ir — T relation) per cent higher than 
the Iron mass accumulated in the ICM through SN activities from 
z = 10. The latter value is, however, well within the scatter in the 
observed distribution. At higher redshift, the observed central value 
in the Tozzi et al. sample of Fe abundance measurements is a factor 
of about 3 higher than the predicted one. 

It is worth noticing that this result does not depend signif- 
icantly upon either the delay time distribution function adopted 
to calculate ri a or the SN CC and la compilations assumed. The 
changes in the total Iron mass accumulated at z = 0.05 are in the 
order of 10 per cent when different cj>(td) are considered and of 
about 40 per cent (i.e. the ratio of the observed and the predicted 



Brief history of the metal accumulation in the ICM 5 



values ranges between 1.0-1.2 and 1.4-1.8; see Table|2j when we 
use the calculations of SN CC explosions in two extreme cases pre- 
sented in Woosley & Weaver (1995) and described at the end of this 
Section. 

Note that, if the SN la rate in number per comoving volume 
and rest-frame year is converted to supernova unit (1 SNu = 1 SN 
per 100 year per 10 10 Lq,b) with an assumed local B-band lumi- 
nosity density of 2 x 10 8 /iLq Mpc -3 which evolves as (1 + z) 1,9 
(see Dahlen et al. 2004 for the caveats in using such conversion 
for rates estimated on cosmological distances), we obtain a local 
(z M 0.05) rate of 0.25 and 0.79 SNu for SN la and CC, respec- 
tively. By integrating the SNu values over the cosmic time, we mea- 
sure a local M Fc /L B of ~ 0.0020 and 0.0010fr? as due to SN 
la and CC, respectively, with a dependence upon the redshift very 
similar to what is shown in Fig.Q The sum of these values is still 
a factor between 3 and 5 lower than the present estimate in galaxy 
clusters: this estimate, however, is affected from the extension of 
the cluster region considered to measure B-band luminosities and 
Iron mass distribution (see discussion on the uncertainties of these 
measurements in De Grandi et al. 2004, Sect. 3. 2). Therefore, even 
though we solve most of the discrepancy between predicted and 
measured Iron mass in typical galaxy clusters through models of 
SN rates that make use of the number of SN per comoving vol- 
ume and rest-frame time, some difficulties remain in recovering 
the observed Mf s /Lb ratio through a direct, but not straightfor- 
ward (given the evolution with redshift of the B-band luminosity 
density), conversion to the number of SN per unit of B-band lu- 
minosity (e.g. Arnaud et al. 1992, Renzini et al. 1993). Regarding 
the latter issue, note that the suggested solution of an increase of 
the SN la rate in E/S0 galaxies by a factor of 5-10 at higher— z is 
partially taken into account in the models presented here (see, e.g., 
Fig.|3j, where the SN la rate is enhanced by a factor of ~ 5 from 
2 = to z ~ 0.8, but then decreases rapidly beyond it. 

We can now extend these considerations to the enrichment his- 
tory of the metals other than Fe. To infer these, we adopt, as done in 
the X-ray analyses presented above, the solar abundances tabulated 
in Anders & Grevesse (1989) and the most recent and widely used 
theoretical metal yields of Nomoto et al. (1997) for (i) SNe CC 
integrated between 10 and 50Mq for a Salpeter IMF (see Naka- 
mura et al. 1999 on the uncertainties related to the adopted mass 
cut that can change the Nickel yield by a factor of 2) and (ii) SNe 
la exploding with fast deflagration according to the W7 model. All 
the numerical values considered here are presented in Table 1. We 
note that the values presented in Nomoto et al. (1997) are plotted 
relative to the solar photospheric abundance in Anders & Grevesse 
(1989). The latter values have been revised to match the meteoritic 
determinations, as summarized in Grevesse & Sauval (1998), and 
require conversion factors of (0.676, 0.794, 1, 1.321, 1) to correct 
the original photospheric abundance for Fe, O, Si, S and Ni, respec- 
tively. We have, however, adopted the Anders & Grevesse values 
for consistency between models and observational constraints. 

Given a metal i, its total mass is estimated as 

Mi = Mi, Ia + Mi,cc 

= 2_j m '- Ia x ria (^) x x dV 

dt,dV 

+ m h cc x rcc(dt) x dt x dV 

dt,dV 

= (M Pe ,z a Y Ia + M Fe ,cc Y cc ) Jf* A : , (4) 



where the synthesized masses for each element of interest, m„ or 
the corresponding abundance yields, Y, are presented in Table 1. 

One direct way to check the consistency of these models is to 
compare the ratio between some of the most prominent elements 
detectable in X-ray spectra and the Iron. To this purpose, we over- 
plot in Fig.|2|the abundance ratios inferred from the models to the 
observational constraints on the O/Fe, Si/Fe, S/Fe and Ni/Fe dis- 
cussed in Section 2. The agreement is remarkably good for what 
concerns the XMM-Newton means and the ASCA Ni/Fe measure- 
ment. A larger departure appears between the model prediction 
and the result for the S/Fe estimate from ASCA. We can evaluate 
the agreement by calculating the x 2 given the observational con- 
straints (with the relative errors) and the predicted values. We ob- 
tain a x 2 = 3.1 (3 degrees of freedom) for the XMM-Newton data, 
which is statistically acceptable [P(x 2 > Xobs) — 0.63] and is, 
however, the largest value measured for any delay time distribution 
functions: the "wide" Gaussian and the e— folding form (see cap- 
tion of Fig.|3j give a x 2 of 2.2 and 2.3, respectively (see Table|2j. 
It is worth noting that the Si/Fe ratio alone contributes about 2.2 
to the total x 2 , owing to the over-predicted amount of Silicon by 
~ 40 per cent. On the other hand, the ASCA measurements provide 
a x 2 of about 10.0 [3 degrees of freedom, P(x 2 > Xobs) = 0.98], 
mainly due to the over-prediction of the S/Fe ratio by 80 per cent. 
We have also considered alternative SN CC and la compilations 
to the Nomoto et al. one, adopted here as reference. Following 
the discussion in Gibson, Loewenstein & Mushotzky (1997) on 
the uncertainties related to the physics of exploding massive stars, 
we consider two extreme conditions for the explosion of SN CC 
as indicated in Woosley and Weaver (1995), namely case A with 
metallicity equals to 10~ 4 solar and case B with solar metallic- 
ity. Both the models provide best-fit results worse than the Nomoto 
et al. one, with a x 2 of 9.3 and 3.1 (case A and B, respectively) 
when compared to the Tamura et al. abundance ratios and 42.5 
and 21.0 when compared to the Baumgartner et al. values (see Ta- 
ble|2j. Moreover, we consider SN la models with a delayed detona- 
tion induced from deflagration at low density layers (model WDD2 
in Nomoto et al. 1997) that seem to reproduce well the observed 
yields in the outer regions of M87 (Finoguenov et al. 2002). This 
model doubles the amount of Si (msi = O.27M0, Y — 1.01) and 
S (ms = 0.17M©, Y = 1.20) and reduces the Ni ejecta by 70 per 
cent (mm = 0.04M Q ,F = 1.40) with respect to W7, amplify- 
ing any mismatch with the observational mean values so that the 
X 2 increases to ~ 8 and 38 when model predictions are compared 
with data from Tamura et al. and Baumgartner et al., respectively. 
We thus exclude the possibility that this model of SN la explosions 
can explain the present results on the overall enrichment history of 
galaxy clusters better than W7. 

We can now trace back the abundance ratios as function of 
redshift (see Fig.[4]and[5j and, given this history of the metals ac- 
cumulation in galaxy clusters, plot in Fig.|6|the relative contribution 
in mass of the SN la products to the total material of each elements 
investigated (i.e. O, Si, S, Fe and Ni) available in the cluster baryon 
budget as function of redshift. 

Locally, 51 per cent of Iron, 1 per cent of Oxygen, 14 per cent 
of Silicon, 21 per cent of Sulfur and 75 per cent of Nichel are pro- 
duced from SN la explosions. These values change only by few 
(< 3) per cent by adopting different delay time distribution func- 
tions and by < 10 per cent when the two other extreme cases of SN 
CC explosion are considered. We conclude that these predictions 
for nearby clusters are particularly stable and robust. As expected 
from the relative increase of SNe CC, the SN la contribution to the 
metal budget decreases significantly at higher redshifts, becoming 
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Table 2. Tests performed on the different SN la (W7, WDD2) and CC (N97, WW95A, WW95B) models and delay time distribution function <f>{t d ) (narrow 
Gaussian, wide Gaussian, exponential form) adopted. The first column indicates the change with respect to the favored model shown in the first row: W7 
deflagration model of SN la products, Nomoto et al. 1997 (N97) compilation of SN CC outputs and a "narrow" Gaussian delay time distribution function. The 
5 tests done are: comparison between the observational constraints (from De Grandi et al. 2004 and Tozzi et al. 2003) on the total amount of Iron mass, Mp c , 
produced at the median redshift of 0.05 and 0.63 in a typical cluster of 5 keV and the predicted value, Mp c , both using a M v j r — T relation from simulations 
(low end) and from observations (high end; 2nd and 3rd column); comparison of the relative change of the Iron abundance with redshift with 49 measurements 
done between redshift 0.3 and 1.3 (Fig.|5] 4th column); comparison of the abundance ratios observed in a sample of nearby clusters with XMM-Newton (O/Fe, 
Si/Fe, S/Fe; Tamura et al. 2004; 5th column) and ASCA (Si/Fe, S/Fe, Ni/Fe; Baumgartner et al. 2005; 6th column). 



models 

(la, CC, <j)(t d )) 


M Fo 
z = 0.05 


z = 0.63 


X 2 (dof) 

vs Fe(z) 


X 2 (dof) 
vsXMM 


X 2 (dof) 
vs ASCA 


(W7, N97, narrow) 


1.29-1.58 


2.58-3.16 


0.64 (5) 


3.1 (3) 


10.0 (3) 


(WDD2, ...) 


1.33-1.62 


2.64-3.22 


0.64 (5) 


8.1 (3) 


37.5 (3) 


(..., WW95A, ...) 


1.43-1.75 


2.94-3.59 


0.66 (5) 


9.3 (3) 


42.5 (3) 


(..., WW95B, ...) 


1.02-1.24 


1.92-2.35 


0.66 (5) 


3.1 (3) 


21.0(3) 


(..., wide) 


1.20-1.47 


2.48-3.03 


0.68 (5) 


2.2 (3) 


8.9 (3) 


(...,..., e-fol) 


1.22-1.49 


2.33-2.85 


0.74 (5) 


2.3 (3) 


9.0 (3) 



responsible for less than half of the metal masses observed locally 
at z ~ 1. More interestingly, the predicted evolution in the Iron 
abundance (Fig.[5J, although significantly steepens up to z ~ 1 
where the expected value is half of the local one, is still consistent 
with observational constraints obtained from a sample of Chandra 
exposures of 49 clusters at redshift > 0.3 (Fig.[5J. 



4 SUMMARY AND DISCUSSION 

We summarize here our main findings on the history of the metals 
accumulation in the ICM. By using the rates of SNe la and CC as 
observed (at z < 1.6 and z < 0.7, respectively) and modelled 
from the cosmic star formation rate derived from UV-luminosity 
densities and IR data sets (Dahlen et al. 2004, Strolger et al. 2004) 
and adopting theoretical yields (as described in Table 1), we infer 
how the metals masses in the ICM are expected to accumulate as a 
function of the redshift. 

We find that these models predict that the total Iron mass ac- 
cumulated in massive galaxy clusters through SN activities is be- 
tween 24 and 75 per cent (235-359 per cent) lower than the Iron 
mass estimated in local (high— z) systems as determined through 
the equivalent width of the emission lines detected in X-ray spectra 
(e.g. De Grandi et al. 2004, Tozzi et al. 2003). This discrepancy is 
reduced by about 20 per cent when the cluster volume used to accu- 
mulate the SN products as a function of time is defined by adopting 
an average virial mass from numerical simulations instead of the 
value extrapolated from the observed X-ray scaling relationships 
(e.g. Ettori et al. 2004b, Arnaud et al. 2005). Because of this kind 
of uncertainty in relating the cluster gas temperature to the asso- 
ciated virial and Iron masses, we conclude that the agreement be- 
tween the expected and measured Iron mass, even though marginal, 
is acceptable within the observed scatter. Moreover, we can repro- 
duce the relative number abundance of the most prominent metals 
detectable in an X-ray spectrum, such as Oxygen, Silicon, Sulfur 
and Nickel with respect to the Iron as estimated for nearby bright 
objects observed both with ASCA (Baumgartner et al. 2005) and 
XMM-Newton (Tamura et al. 2004). We also show that the predicted 



evolution in the Iron abundance, which should decrease by a factor 
of 2 at z k, 1 with respect to the local value, is in good agreement 
with the current observational constraints obtained from a sample 
of 49 clusters at z > 0.3 (see Fig.[5). 

By using these models to describe the ICM enrichment, we 
can infer the relative SN contribution to the amount of elements 
present in galaxy clusters as function of the cosmic time. At in- 
creasing redshift, the products from SNe CC become dominant, 
owing to the steep rise of their relative rate with respect to SNe 
la. The transition occurs between z = 0.5 and 1.4, with an en- 
hancement of the a— elements with respect to Fe (e.g. O/Fe and 
Si/Fe ratios increase by a factor of more than 2) and a drastic de- 
crease of the Ni/Fe ratio (see Fig.|4j. Then, the fractions of metals 
mass due to SN la outputs, which are locally about 51, 75, 1, 14 
and 21 per cent of the total for Iron, Nickel, Oxygen, Silicon and 
Sulfur, respectively, halve at z ~ 1 (Fig. [6), almost independently 
from the adopted delay time distribution functions and models for 
SN CC explosions. When the assumed SN rates (number per co- 
moving volume per rest-frame year) are converted to units of the 
B-band luminosity (but see caveats in Dahlen et al. 2004), which 
is well mapped only in nearby galaxy clusters, we show that the 
expected total Mf c /Lb is still below the local measurement by a 
factor between 3 and 5 and that SN la metal production contributes 
by ~70 per cent. We conclude that this well-known (e.g. Arnaud et 
al. 1992, Renzini et al. 1993) underestimate of the local M^/Lb 
value cannot be explained with an increase in the SN rates at higher 
redshift, according to the present models. 

The changes by a factor of 2 in the abundance ratios at z > 0.5 
arising from the predominance of the enrichment through SNe CC 
might be investigated in the near future with X-ray spectroscopi- 
cally resolved metal abundance estimates in high-redshift galaxy 
clusters. We are attempting this with our Chandra and XMM- 
Newton data set of high— z objects (Tozzi et al. 2003, Ettori et 
al. 2004a), with expected relative uncertainties on the abundance 
ratios of larger than ~ 20 per cent at the la level (see, for ex- 
ample, Fig. |5j- These X-ray satellites offer the best compromise 
available at present between field-of-view, effective area, and the 
spatial and spectral resolution required to pursue such a study. Only 
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Figure 5. (Upper panel) Distribution of the metal abundance as a func- 
tion of redshift for a sample of Chandra exposures of 49 galaxy clusters at 
z > 0.3 and gas temperatures between 3.4 and 15.5 keV (Balestra, Tozzi, 
Rosati, Ettori et al. in preparation, which will present the metallicity mea- 
surements for a sample larger by a factor of ~ 3 than the one presented 
in Tozzi et al. 2003, with upgraded calibration files and extended statistical 
analysis). The diamonds refer to the median distribution of 10 Iron abun- 
dance measurements per redshift bin. (Lower panel) Predicted and observed 
Iron abundance as a function of redshift. The diamonds refer to the median 
distribution shown in the upper panel. The normalization of the predicted 
evolution is here fixed equal to the first measured bin at Z Ri 0.4. The lines 
refer to different SN explosion models and delay time distribution functions 
with respect to the reference values indicated by the solid line (dashed line: 
Woosley & Weaver SN CC models, case A -thin line- and B -thick line-; 
dotted line, from the thinnest to the thickest line: WDD2 SN la model, wide 
Gaussian and e— folding form). 




redshift 

Figure 6. Metal mass fraction resulting from SN la ejecta as a function of 
redshift. 

with XEUS 1 and its sensitivity greater by two orders of magnitude 
than that of XMM-Newton and a spectral resolution of the order of 
10 eV or less, will it be possible to investigate with significantly 
higher accuracy the metal budget of the ICM in high— z systems. 
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